function preset()
{
    #=================================
    fig_name=bathy_A3
    fig_fmt=png
    #=================================
    #TOC=../../gmt/TOC_rupture/toc_area.grd
    Locking=../inteface/Schurr2014/Schurr_2014_Nature_locking.grd
    main=../interface/Schurr2014/Schurr_2014_Nature_main.grd
    bathy=../../gmt/Geersen_2018_Tectonics.nc
    after=../interface/Schurr2014/Schurr_2014_Nature_after.grd
    aftershock=../interface/flo_aftershock/after_latlongdep.txt
    bathy_merge=../../gmt/P1_3D_merge.grd
    
    # 透明度
    alpha_profile=50
    #图片大小
    fig_x=10
    fig_y=8
    fig_z=2
    #数据范围
    zmin=-20000
    zmax=1000
    lon_min=-72
    lon_max=-69.5
    lat_min=-22.2
    lat_max=-18.4
    #range
    R=$lon_min/$lon_max/$lat_min/$lat_max
    J=m7.8c #1°对应 1cm,那么如果要满足 横向 16c 的话 1°要对应 6.4cm
    #colorbar
    gmt makecpt -Crainbow -T0/5 > rupture.cpt
    gmt makecpt -Ccyclic -T7000/38000 > toc_area.cpt
    oceanic=../../gmt/P1_basemap_shade/oceanic.cpt
    #intens=../bathy_A3/intens.grd
}

preset

gmt begin $fig_name $fig_fmt E2000
    #gmt grdcut $bathy -R$R -V -Ggrdcut.grd
    #gmt grdgradient grdcut.grd -Nt1 -A30 -Gintens.grd -V
    gmt basemap -R$R -J$J -Ba0.5f0.5 -BWSne
    gmt grdimage $bathy -R$R -J$J -C$oceanic -Iintens.grd 
    gmt psscale -DjTL+w2.2c/0.2c+o1.3i/0.2i -C$oceanic -Ba4000f2000 # -Dx3.5c/13c+w2.2c/0.2c
    gmt coast -R$R -J$J -W0.25p,black -Ggrey
    awk '$3 ~ /Arica/ {print $1,$2}' ../P1_fig/Chile_cities.txt > tmp.txt
    gmt plot tmp.txt -R$R -J$J -W0.1c,black -Ss-0.1c
    awk '$3 ~ /Iquique/ {print $1,$2}' ../P1_fig/Chile_cities.txt > tmp.txt
    gmt plot tmp.txt -R$R -J$J -W0.1c,black -Ss-0.1c
    awk '$3 ~ /Tocopilla/ {print $1,$2}' ../P1_fig/Chile_cities.txt > tmp.txt
    gmt plot tmp.txt -R$R -J$J -W0.1c,black -Ss-0.1c
    for v in 04 23 25 17 32 11 15 16A 30 13 #dipping line and strike line
    do
        gmt psxy ../../gmt/P1_basemap_shade/MGL1610MC${v}_new.input -R$R -J$J -W0.06c,black 
        gmt pstext ../../gmt/P1_basemap_shade/MGL1610MC${v}.annotation -R$R -J$J -F+a-90+f5p,Helvetica-Narrow-Bold -D0c/-0.3c 
        gmt psxy ../../gmt/P1_basemap_shade/MGL1610MC${v}.coord -R$R -J$J -W0.03c,red -Sc-0.05c 
    done
    gmt psbasemap -LjRT+o1c/0.5c+c-71/-21+w30k+f --MAP_SCALE_HEIGHT=0.3c
    gmt inset begin -DjTL+w1i+o0.15i/0.1i 
        gmt pscoast -Rg -JA-70.2/-20/2.5c -Bg -Dl -A1000 -Ggrey -Swhite 
    gmt inset end
gmt end show